Fermi-Bose transformation for the time-dependent Lieb-Liniger gas 
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Exact solutions of the Schrodinger equation describing a freely expanding Lieb-Liniger (LL) gas 
of delta-interacting bosons in one spatial dimension are constructed. The many-body wave function 
is obtained by transforming a fully antisymmetric (fermionic) time-dependent wave function which 
obeys the Schrodinger equation for a free gas. This transformation employs a differential Fermi-Bose 
mapping operator which depends on the strength of the interaction and the number of particles. 
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Noncquilibrium phenomena in quantum many-body 
systems are among the most fundamental and intrigu- 
ing phenomena in physics. One-dimensional (ID) inter- 
acting Bose gases provide a unique opportunity to study 
such phenomena. In some cases, the models describing 
these systems [H, 0, H| allow to determine exact time- 
dependent solutions of the Schrodinger equation 0, Q 
providing insight beyond various approximations, which 
is particularly important in strongly correlated regimes. 
These ID systems are experimentally realized with atoms 
tightly confined in effectively ID waveguides [1, 0, [1], 
where nonequilibrium dynamics is considerably affected 
by the kinematic restrictions of the geo metry Q, while 
quantum effects are enhanced (jl. HoL Today, exper- 
iments have the possibility to explore ID Bose gases for 
various interaction strengths, from the Lieb-Liniger (LL) 
gas with finite coupling [fj, |8(] up to the so-called Tonks- 
Girardeau (TG) regime of "impenetrable-core" bosons 
0, HI- However, most theoretical studies of the ex- 
act time-dependence address the TG regime (see, e.g., 

Refs. HEtHQHESHEl)- Inthis limit > the com - 

plex many-body problem is considerably simplified due to 
the Fermi-Bose mapping property [4j where dynamics fol- 
lows a set of uncoupled single-particle (SP) Schrodinger 
equations (4j. It is therefore desirable to employ an effi- 
cient method for calculating the time-evolution of a LL 
gas with finite interaction strength. 

In 1963, Lieb and Liniger [l[ presented, on the ba- 
sis of the Bcthe ansatz, a solution for a homogeneous 
Bose gas with (repulsive) S- function interactions, for ar- 
bitrary interaction strength c; periodic boundary con- 
ditions were imposed. This system was analyzed by 
McGuire on an infinite line with attractive interactions 
3]. The renewed interest in ID Bose gases stimulated 
recent studies of static LL wave functions OHHH in- 
cluding a LL gas in box confinement 2l|. Besides the 
wave functions, the correlations of a LL system with fi- 
nite coupling [2l, [H, 0, M, El, M, El, M, H [HI pro- 
vide a link to many observables and were analyzed by 
using various techniques, including the inverse scattering 
method [3, 25 , HO, Ei , 1/c expansions [3| relying on 



the analytic results in the TG regime 1 32] , and numerical 
Quantum Monte Carlo techniques [28(. Regarding dy- 
namics, a full numerical study of the irregular dynamics 
in a mesoscopic LL system was presented in [33| . In Ref . 
[H, Girardeau has shown that phase imprinting by light 
pulses conserves the so-called cusp condition imposed by 
the interactions on the LL wave functions, and suggested 
to use time-evolving SP wave functions to analyze the 
subsequent dynamics. However, as pointed out in Ref. 
H, the presented scheme does not obey the cusp condi- 
tion during the evolution which limits its validity. This 
situation can be remedied by using an ansatz which obeys 
the cusp condition at all times by construction [2~il l34j. 

Here we construct exact solutions for the freely ex- 
panding LL gas with localized initial density distribution. 
This can be achieved by differentiating a fully antisym- 
metric (fermionic) time-dependent wave function, which 
obeys the Schrodinger equation for a free Fermi gas 34j ; 
the employed differential operator depends on the inter- 
action strength c and the number of particles. When c — > 
00, the scheme reduces to Girardeau's time-dependent 
Fermi-Bose mapping valid for "impenetrable-core" 
bosons. 

We consider the dynamics of N indistinguishable 6- 
interacting bosons in a ID geometry [l| . The Schrodinger 
equation for this system is 
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N 

■E 

i=l 



2c5{xi 

l<i<j<N 



Xj)lp B , (1) 



where i/)b{%i, ■ ■ ■ , %N, t) is the many-body wave function, 
and c quantifies the strength of the interaction (for con- 
nection to physical units see, e.g., [B|). The a;-space is in- 
finite (we do not impose any boundary conditions) , which 
corresponds to a number of interesting experimental situ- 
ations where the gas is initially localized within a certain 
region of space and then allowed to freely evolve. This is 
relevant for free expansion [l3|> [IJ, 2 IB] or interference 
of two initially separated clouds during such expansion 
[13] , etc. Due to the Bose symmetry, it is sufficient to ex- 
press the wave function ips in a single permutation sec- 
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tor of the configuration space, R\ : x\ < X2 < ■ ■ ■ < xn- 
Within ipB obeys 



N 



where the primed operator CK - +1 = 0/Bj.j + i omits the 
factor Bjj + i as compared to O. The auxiliary function 
can be written as 



idxij B /dt = -y d 2 ip B /dxt, 



i=l 



(2) 
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while interactions impose boundary conditions at the 
borders of Ri []]]: 
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This constraint creates a cusp in the many-body wave 
function when two particles touch, which should be 
present at any time during the dynamics. 

In the TG limit (i.e., when c — > oo) the cusp condi- 
tion is iPb(xi, ■ ■ .,Xj,Xj+i, . . .,x N ,t)\ Xj+1=Xj = @, 0], 
which is trivially satisfied by an antisymmetric fermionic 
wave function ippixi, . . . , xn, i); thus ipB = V'F within 
which is the famous Fermi-Bose mapping @,5]. In 
many physically interesting cases, tpF can be constructed 
as a Slater determinant 

i/> F (xi t ...,x N ,t) = {N\)-idct^ m (x 3 ,t)}^ 3=1 . (4) 

Since ipB = ipF within R±, ipp must obey idipF/dt = 
— d 2 tpF / dx 2 p which implies that the (orthonormal) 

SP wave functions <j) m (xj , t) evolve according to 

id(f> m /dt = ~d 2 <j) m /dx 2 - (5) 

m = 1, ...,N. Thus, in the TG limit, the complexity 
of the many-body dynamics is reduced to solving a sim- 
ple set of uncoupled SP equations, while the interaction 
constraint ([3]) is satisfied by the Fermi-Bose construction. 

The simplicity and success of this idea motivates us to 
choose an ansatz which automatically satisfies constraint 
([3]) for any finite c 24|, |34( . For this, define a differential 
operator 
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Bij - 



1 



d 

dXj 



'3i 



JL 

dx,. 



It is straightforward to verify that the operator 
Bj + i jBjj + iO'j - +1 in front of tpp is invariant under 
^p B = o. (3) the exchange of Xj and Xj+i [24]. On the other hand, 
the fermionic wave function ipp is antisymmetric with 
respect to the interchange of Xj and Xj+i. Thus, 
^AUxC^ij ■ ■ ■ , x ji x j+ii ■ ■ ■ i x n, t) is antisymmetric with 
resp ect to the interchange of Xj and Xj+i, which leads to 
HS] 

1pAUx(xi, ■ ■ ■ ,Xj,Xj + i, . . . ,x N ,t)\ Xj+1=Xj = 0. (11) 

This is fully equivalent to the cusp condition ([3]), 
Bj+ijipB\x j+ i=Xj = [ISjS]- Thus, the wave function 
dSJ obeys constraint ((3j) by construction. 

In order to exactly describe the dynamics of LL gases, 
the wave function (JSj> should also obey Eq. ([2|) inside R\ . 
From the commutators [d 2 /dx 2 , O] = and [id/dt, O] = 
it follows that if ipF is given by Eq. ([3]) and the 4> m (xj , t) 
obey Eq. (J5]), then ipg obeys Eq. (J5J. Note that for c — » 

00, one recovers Girardeau's Fermi-Bose mapping 0,0], 

1. e., 6 = 1. 

Let us utilize this formalism to describe the dynamics 
of a freely expanding LL gas. Suppose that for t < 
the system is confined by an external potential V(x) 
and is in its ground state, before at t = the poten- 
tial is suddenly switched off. In order to find the ex- 
act form of the initial condition, we have to solve the 
static Schrodinger equation for the LL gas in the poten- 
tial V(x). By using the above formalism, we express the 
initial state as ^bo = NcOipFo, which should (within R\) 
obey Y,j HjipBQ = E B o^bo or, equivalently, 

Oj^HjipFo - [O^Hj^FO = OEbqiPfo- (12) 



(6) 



It can be shown that the wave function 
%j)B = A/" c O-0f (inside Ri), 



(7) 



(8) 



where j\f c is a normalization constant, obeys the cusp 
condition ((3]) by construction 24, 34|: Consider an aux- 
iliary wave function 



4>AUX(X1, ■ ■ - ,X N ,t) 



B 



B j+ ijB jtj+1 d'^F, (9) 



Here, Ebo is the ground-state energy, and Hj = 
—d 2 /dx 2 + V(xj) is the SP hamiltonian. Eq. (fl"2"j) shows 
that, due to the nonvanishing commutator [0,^2j Hj] = 

[O, J2j V{xj)], operating with O on the fermionic ground 
state in the trap does not give the bosonic ground 
state. However, for sufficiently strong interactions 
and/or weak and slowly varying potentials, we can ap- 
proximate [0,Y2j^i x j)] ~ 0- I n t nc TG limit, the 
commutator vanishes identically. Thus, for sufficiently 
strong interactions, the ground state is approximated by 
-050 =N c 6det[(f> m (xj,Q)]^ j=1 /y/M, where </> m (xj,Q) is 
the mth eigenstate of the SP hamiltonian. 

In what follows we study free expansion from such 
an initial condition, which describes a LL gas with a 
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localized density distribution. Even though ipso can 
not be interpreted as a ground state for weak inter- 
actions, the free expansion from 7p B o is calculated ex- 
actly for all values of c. We illustrate the properties of 
ipBo = Af c O det[4> m (xj , 0)]^ j =1 /\/N\ for the harmonic 
potential V(x) = v 2 x 2 /4, with v = 2. Fig. CD (left col- 
umn) displays the section |V>so(0, Xi, X3) | 2 of the prob- 
ability density, for N = 3 particles and three values of 
c. We clearly see that, as the interaction strength in- 
creases, the initial state becomes more correlated. Given 
that one particle is located at zero, for c = 1 there is a 
considerable probability that the other two particles are 
to the left or to the right of the first one, i.e., their posi- 
tions are weakly correlated with that of the first particle. 
However, for larger c, if one particle is at zero, it is more 
likely that the other two particles are on opposite sides 
of the first one, and their distance grows with increasing 
interaction strength. 

When the harmonic potential is turned off, the evolu- 
tion of the SP states <j) m {x,t) is known exactly (see, e.g., 
Ref. [H): 4> m (x,t) = cf> m (x/b(t),0) exp[ix 2 &' (t)/(46(t))- 
iE m T(t)]/ 'y/b(t), where E m is the energy of the mth 
SP eigenstate <j> m (x,0), b(t) = \/l + t 2 u 2 , and r(t) — 
arctan(z4)/V. We can make use of the expression for the 
ground state of a TG gas in harmonic confinement [35} 
to calculate Vfo- Employing Eq. the evolution of 
the many-body wave function t/jg can then be formally 
expressed (within Ri) as 



Oe 



ip B =M(c,v,N)b{t)-- 
T^£f = i[^7&(*)] 2 TJ 

l<i<j<N 



i X 3 



~r(t) 

Xi), 



(13) 



where 7V(c, N) is a normalization constant, evaluat- 
ing to J\f — 1 for c — > 00. The action of the operator 
O yields lengthy expressions already for a few particles, 
and particular examples will be given elsewhere. The 
asymptotic form of Eq. (|13| is given by limt_»oo ips /ipF ot 

{on Ki<j<N( X J x i)}/ T\l<i<j<N( X j Xi )- Although 

Eq. (|13p provides an exact wave function for the time- 
dependent LL gas, it is desirable to calculate the evo- 
lution of observables such as the SP density p(x, t) = 
N J dx2 ■ ■ ■ (Ixn\iPb(x, X2, ■ ■ ■ , Xn, t)\ 2 - This task is com- 
plicated by the many-fold integral. However, we can find 
the evolution of p(x, t) numerically for small numbers of 
particles. Fig. [T] (right column) displays the evolution of 
the SP density for three different values of c. For larger 
c, the initial SP density exhibits typical TG-fermionic 
properties, characterized by N small separated humps 
[J, |l2j . For all values of c, the SP density acquires such 
humps during free expansion indicating that the system 
becomes more correlated in time, which is in accord with 
the fact that the LL gas becomes strongly correlated for 
lower densities 

Even though the employed approach is valid at any 
interaction strength, it is particularly useful for the 
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FIG. 1: (color online) The correlation properties of the initial 
state and the evolution of the SP density for various interac- 
tion strengths, for N = 3 particles. (Left column) The proba- 
bility density \4>bo{0, %2, %3)\ 2 for c = 1,3, 10. (Right column) 
The SP density at t = (black solid line), t = 0.5 (red dotted 
line), and t = 1 (blue dot-dashed line), for c = 5, 10, 50. 



strongly interacting gas: First, the operator O can be 
hierarchically organized into orders 1/c , O = 1 + 

J2k=i ° k Ok- By keeping the terms of order 1/c, 
we obtain the first-order correction to the TG gas. In 
this approximation, the form of the operator is consider- 
ably simplified, O^l + r'Oi, where 0\ = ^2k = i(2k — 
TV — l)d/dxk- The wave function reads, within 



1pB=1pF + 



N 

E 



2k- N 



-Act[A k n 



N 

mj lm,j- 



(14) 



where A k nj = 4> m (xj,t) for j ^ k, and A k mk = 
d<t> m {xk-~t) I d x k- The numerical calculation of the wave 
function (I14|) is not a difficult task even for a fairly 
large numbers of particles. Second, in this regime, 
[0,J2jHj] = [0,Y^,jV(Xj)] si is a reasonable ap- 
proximation. For example, with O ~ 1 + c 0\ and 
V(x) = i/ 2 x 2 /4, [0,J2j v ( x j)} is 01 order i/ 2 /4c, i.e., for 
v 2 < 1/c the commutator is of order 1/c 2 or less. Thus, 
the approach can be used to characterize time-dependent 
and static LL gases in various trapping potentials in the 
strongly correlated regime, but below the TG gas limit. 

For completeness, let us briefly discuss time-evolving 
states ipB with periodic boundary conditions as in Ref. 
l]. Any time-evolving state ipB can be written as a 
superposition of eigenstates iphh,c(xi, . . . ,xn), where £ 
denotes all quantum numbers necessary to describe one 
eigenstate. LL eigenstates can be written as = 



4 



OA/" e det[e ife ™^']^ J=1 [24j,|34{. If periodic boundary con- 
ditions are imposed as in Ref. [lj, the quasimomenta 
kj must obey a set of coupled transcendental equations 
and depend on c [H, Time-evolving states ipB can be 
written as a superposition of LL eigenstates 



-iEct 



where the coefficients &(£) are fixed by the initial condi- 
tions. These coefficient are in practice hard to calculate 
given the initial state due to the many-fold integrations 
that need to be performed. 

It should be emphasized that the Fermi-Bose trans- 
formation employed here differs from the fermion-boson 
duality discussed by Cheon and Shigehara [36| (see also 
[37l] ) . because it transforms a noninteracting fermionic 
wave function into a wave function describing LL gas. 
Using Ref. [36[ it can be shown that the approach used 
here can also be applied to construct wave functions for 
a time-dependent Fermi gas with finite-strength interac- 
tions. 

In conclusion, we have constructed exact solutions for 
the freely expanding LL gas with localized initial density 
distribution. Wave functions are obtained by differenti- 
ating a fully antisymmetric (fermionic) time-dependent 
wave function, which obeys the Schrodingcr equation for 
a free Fermi gas. For a number of physically interesting 
situations (e.g., free expansion), by using the operator 
O, the state of the system can be derived from N SP 
time-dependent states, as anticipated in Ref. [B|. The 
construction of LL wave functions for various external 
potentials V(x), and the derivation of correlation func- 
tions within the employed formalism is the subject of 
ongoing work. 
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